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Using the Hellmann-Feynman theorem we analyze the contribution of the different terms of the 
nucleon-nucleon interaction to the nuclear symmetry energy Esym and the slope parameter L. The 
analysis is performed within the microscopic Brueckner-Hartree-Fock approach using the Argonne 
V18 potential plus the Urbana IX three-body force. We find that the main contribution to Esym 
and L is due to the tensor component of the nuclear force. 
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The nuclear symmetry energy, defined as the differ- 
ence between the energies of neutron and symmetric mat- 
ter, and in particular its density dependence, is a cru- 
cial ingredient to understand many important properties 
of isospin-rich nuclei and neutron stars [l|-|3|. Exper- 
imental information on the density dependence of the 
symmetry energy Esym{p) below, close to and above 
saturation density pq can be obtained from the analy- 
sis of data of isospin difusion measurements j^], giant 
and pygmy resonances Q, isobaric analog states Q, 
isoscaling i8| or meson production in heavy ion colli- 
sions p, 10]. Accurate measurements of the neutron skin 

thickness SR = \/{Rn) ^ \ {Rp) in heavy nuclei, via 
parity- violating electron scattering experiment s fill [l2| 
or by means of antiprotonic atom data [IJ, [ij, can 
also help to constraint Esym{p]i, since its derivative its 
strongly correlated with SR [la]- Additional informa- 
tion on Esym {p) can be extracted from the astrophysical 
observations of compact objects which open a window 
into both the bulk and microscopic properties of nuclear 
matter at extreme isospin asymmetries [3|. In particu- 
lar, the characterization of the core-crust transition in 
neutron stars [16l - l2g |. or the analysis of power-law cor- 
relations, such as the relation between the radius of a 
neutron star and the equation of state [2l| can put strin- 
gent constraints on Egymip)- Theoretically Esym{p) has 
been determined using both phenomenological and mi- 
croscopic many-body approaches. Phenomenological ap- 
proaches, either relativistic or non-relativistic, are based 
on effective interactions that are frequently built to re- 
produce the properties of nuclei [22| . Since many of such 
interactions are built to describe systems close to the 
symmetric case, predictions at high asymmetries should 
be taken with care. Skyrme-Hartree-Fock [23] and rela- 
tivistic mean field [24] calculations are the most popular 
ones among them. Microscopic approaches start from 
realistic nucleon-nucleon (NN) interactions that repro- 
duce the scattering and bound state properties of the 
free two-nucleon system and include naturally the isospin 
dependence [25|. The in- medium correlations are then 
built using many-body techniques that microscopically 



account for isospin asymmetric effects such as, for in- 
stance, the difference in the Pauli blocking factors of 
neutrons and protons in asymmetric matter. Among 
this type of approaches the most popular ones are the 
Brucckncr-Bcthe-Goldstone fBB Gl ^q and the Dirac- 
Brueckner-Hartree-Fock (DBHF) [23| theories, the vari- 
ational method [28[, the correlated basis function (CBF) 
formalism [2^, the self-consistent Green's function tech- 
nique (SCGF) [30] or, recently, the Yiowk approach [3l| . 
Nevertheless, in spite of the experimental [33 and the- 
oretical [33] efforts carried out to study the properties 
of isospin-asymmetric nuclear systems, Esymip) is still 
uncertain. Its value Egym at saturation is more or less 
well established {^ 30 MeV), and its behavior below sat- 
uration is now much better known [34| . However, for 
densities above po, Esymip) is not well constrained yet, 
and the predictions from different approaches strongly 
diverge. Why Esym{p) is so uncertain is still an open 
question whose answer is related to our limited knowl- 
edge of the nuclear force, and in particular of its spin 
and isospin dependence |35l - l43 ]. 

In this letter we analyze the contribution of the differ- 
ent terms of the NN interaction to Egym and the slope pa- 
rameter L = 3 po {dEsymip)/dp)pg. The analysis is car- 
ried out with the help of the Hellmann-Feynman theorem 
[431 ] within the framework of the microscopic Brueckner- 
Hartree-Fock (BHF) approa ch [261] . We employ the Ar- 
gonne V18 (Avl8) potential 1441 supplemented with the 
Urbana IX three-body force [4^ which for the use in 
the BHF approach is reduced to an effective two-body 
densit y-d ependent force by averaging over the third nu- 
cleon [4y]. We find that the tensor term of the nuclear 
force gives the largest contribution to both Esym and L. 

The BHF approch is the lowest order of the BBG 
many-body theory [26j . In this theory, the ground state 
energy of nuclear matter is evaluated in terms of the so- 
called hole-line expansion, where the perturbative dia- 
grams are grouped according to the number of indepen- 
dent hole- lines. The expansion is derived by means of 
the in-medium two-body scattering G-matrix. The G- 
matrix, that takes into account the effect of the Pauli 
principle on the scattered particles and the in-medium 
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TABLE I; Kinetic (T) and potential {V) contributions to 
Emm, Esm, Esym and L. Units are given in MeV. 



potential felt by each nucleon, has a regular behavior 
even for strong short-range repulsions, and it describes 
the effective interaction between two nucleons in the pres- 
ence of a surrounding medium. In the BHF approach, the 
energy is given by the sum of only two-hole-line diagrams 
including the effect of two-body correlations through the 
G-matrix. It has been shown by Song et al., I47| that the 
contribution to the energy from three-hole-line diagrams 
(that account for the effect of three-body correlations) is 
minimized when the so-called continous prescription [48| 
is adopted for the in-medium potential, which is a strong 
indication of the convergence of the hole-line expansion. 
We adopt this prescription in our calculation. 



The BHF approach does not give access to the separate 
contributions of the kinetic and potential energy in the 
correlated many-body state, because it does not provide 
the correlated many-body wave function |^). However, it 
has been recently shown [49[ that the Hellmann-Feyman 
theorem [43] can be used to calculate the ground state ex- 
pectation values of both contributions from the derivative 
of the total energy with respect to a properly introduced 
parameter. Writing the nuclear matter Hamiltonian as 
H = T -\- V , and defining a A dependent Hamiltonian 
H{X) ~ T -\- \V, the expectation value of the potential 
energy is given as 
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Then, the kinetic energy contribution (T) can be ob- 
tained simply by subtracting {V) from the total energy 
E. 

In Table U we show the kinetic and potential contribu- 
tions to the energy of neutron matter Ej^m, symmetric 
matter Esm, Esym and L at saturation (po — 0.187 fm^^ 
in our calculation) . We note that the kinetic contribution 
to Esym is very small and negative. This is in contrast 
with the result for a free Fermi gas (FFG), whose con- 
tribution at pq is ~ 14.4 MeV. A similar result has been 
recently found by Xu and Li [4l| . According to these au- 



thors, this is due to the strong isospin-dependence of the 
short-range NN correlations (SRC) induced by the tensor 
force. They have shown, in fact, that the increase of the 
kinetic energy of symmetric matter due to SRC is much 
larger than that of neutron matter, the kinetic part of 
the symmetry energy becoming then negative. We also 
note that the kinetic contribution to L is smaller than 
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TABLE II: Partial wave decomposition of the potential part 
of Enm, Esm, Esym and L. Units are given in MeV. 



the corresponding one of the FFG (L^^<^ - 29.2 MeV). 
The major contribution to both Esym and L is due to 
the potential part. Note that, in fact, this contribution 
is practically equal to the total value of Esym and it rep- 
resents '- 78% of L. 



Tables |TT] and IIIII show the partial wave, and the spin 
(S) and isospin (T) channel decompositions of the poten- 
tial part of Enm, Esm, Esym and L at po- Contributions 
up to J = 8 have been considered. We observe that the 
spin-triplet {S = 1) and isospin-singlet (T = 0) channel, 
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TABLE III: Spin (S) and isospin (T) channel decomposition 
of the potential part of Emm, Esm, Esym and L. Units are 
given in MeV. 



and in particular the ■^S'l wave, gives the largest contri- 
bution to both Egyjn and L. This is due, as we explicitly 
show in the following, to the effect of the tensor compo- 
nent of nuclear force that dominates the potential contri- 
bution to the symmetry energy and L, mainly through 
the ^5*1 —^ Di channel. Note that this channel, which 
gives the major contribution to the energy of symmetric 
matter, does not contribute to neutron matter. Note also 
that isospin-triplet (T = 1) channels give similar contri- 
butions to both Enm and Esym which almost cancel out 
in Esym- Similar arguments have been pointed out by 
other authors 1351-1421. 



Next, we analyze the role played by the different terms 
of the nuclear force, particularly the one of the tensor, 
in the determination of Esym and L. To such end we 
apply the Hellmann-Feynman theorem to the separate 
components of the Avl8 potential and the Urbana IX 
three-body force. The Avl8 potential has 18 components 
of the form Vp{rij)Ofj with 
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TABLE IV: Separate contributions to Enm, Esm, Eaym and 
L from the various components of the Avl8 potential (de- 
noted as {Vi)) and the reduced Urbana force (denoted as (Ui)). 
Units are given in MeV. 
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Oij , Jij {Ti ■ Tj), L ■ O, L ■ b{Ti ■ Tj), L , 

i^(Tj • Tj),L'^{di ■ dj), L^{di ■ dj){n ■ fj), {l ■ sY, 

{L ■ S) {n ■ f,j), Ti-j , {di ■ dj)Tij , Sij Tij , {Tzi + Tz- ) 

being Sij the usual tensor operator, L the relative orbital 
angular momentum, S the total spin of the nucleon pair, 
and Tij = Sr^iTzj — t.; • r, the isotensor operator defined 
analogously to Sij. Note that the last four operators 
break the charge independence of the nuclear interaction. 



As we said above, the Urbana IX three-body force is 
reduced to an effective density-dependent two-body force 
when used in the BHF approach. For simplicity, in the 
following we refer to it as reduced Urbana force. This 
force is made of 3 components of the type Up{rij , p)0^j 

where 0'^~ ' = l,{di-dj){fi-fj), Sij (f, -fj), introducing 
additional central, ar and tensor terms (see e.g., Baldo 
and Ferreira in Ref. [4y| for details). 

The separate contributions to Ejsim, Esm, Esym and L 
from the various components of the Avl8 potential and 
the reduced Urbana force are given in Table IIVI The 
contribution from the tensor component to Esym and L 



(contributions (Vsij) and (V5. .(?..?.)) from the Avl8 po- 
tential, and {Us..tf..fA) from the reduced Urbana force) 
is 36.056 MeV and 69.968 MeV, respectively. These re- 
sults clearly confirm that the tensor force gives the largest 
contribution to both Esym and L. The contributions 
from the other components are either negligible, as for in- 
stance the contribution from the charge symmetry break- 
ing terms {{Vt,^), (V(5=..ff^)T„), {Vs,,t,,) and (V(^„+r,p)), 
or almost cancel out. 

In summary, using the Hellmann-Feynman theorem we 
have evaluated the separate contribution of the different 
terms of the nuclear force to the nuclear symmetry energy 
Esym and the slope parameter L. Our study has been 
done within the framework of the BHF approach using 
the Avl8 potential plus an effective density-dependent 
two-body force deduced from the Urbana IX three-body 
one. Our results show that the potential part of the nu- 
clear Hamiltonian gives the main contribution to both 
Esym and L. 



The kinetic contribution to Esym is very 



small and negative in agreement with the recent results of 
Xu and Li [41| . We have performed a partial wave, and a 
spin-isospin channel decomposition of the potential part 
of Esym and L, showing that the major contribution to 
them is given by the spin-triplet {S = 1) and isospin- 
singlet (T = 0) channel. This is due, as we have ex- 



plicitly shown, to the dominant effect of the tensor force 
which gives the largest contribution to both Egy^ and L. 
In conclusion, our results confirm the critical role of the 
tensor force in the determination of the symmetry energy 
and its density dependence. 
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